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Abstract. We present the results of a study aimed at de- 
termining the 12 C/ 13 C ratio in two samples of planetary 
nebulae (PNc) by means of millimeter wave observations 
of 12 CO and 13 CO. The first group includes six PNe which 
have been observed in the 3 He + hyperfine transition by 
Balser et al. (1997); the other group consists of 22 neb- 
ulae with rich molecular envelopes. We have determined 
the carbon isotopic ratio in 14 objects, 9 of which are 
new detections. The results indicate a range of values of 
12 C/ 13 C between 9 and 23. We estimate the mass of the 
progenitors of the PNe of our sample and combine this in- 
formation with the derived 12 C/ 13 C isotopic ratios to test 
the predictions of stellar nucleosynthesis models. 

We find that the majority of PNe have isotopic ratios 
below the values expected from current standard asymp- 
totic giant branch models in the mass range of interest. We 
suggest that the progenitors of the PNe must have un- 
dergone a non-standard mixing process during their red 
giant phase and/or asymptotic giant phase, resulting in 
a significant enhancement of the 13 C abundance in the 
surface layers. Our study confirms a similar behaviour in- 
ferred from spectroscopic observations of field population 
II stars and globular cluster giants, and extends it to the 
final stages of stellar evolution. Finally, we discuss the im- 
plications of our results on models of Galactic chemical 
evolution of 3 He and 12 C/ 13 C. 
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1. Introduction 

In the PN phase, stars more massive than 0.8 Mq re- 
turn to the ISM material that has been processed in the 
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stellar interior. This matter mixes with the surrounding 
medium and modifies the original abundances of elements. 
The contribution of PNe to the Galactic chemical evo- 
lution is particularly important for 3 He which, together 
with deuterium, plays a fundamental role in testing the 
standard Big Bang nucleosynthesis model. While the evo- 
lution of deuterium is well understood, that of 3 He still 
encounters serious problems which cast doubts on the use- 
fulness of this isotope as a test of Big Bang nucleosyn- 
thesis models (e.g. Galli et al. 1995; for a different opin- 
ion on deuterium see Mullan &; Linsky 1999). Observa- 
tions of 3 He toward PNe and Hn regions have resulted 
in abundances that differ by almost two orders of mag- 
nitude: 3 He/H - 1(T 3 in PNc (Rood et al. 1992, Balser 
et al. 1997), and 3 He/H - 1CT 5 in Hn regions (Balser et 
al. 1994, Rood et al. 1995). The latter value is also rep- 
resentative of the 3 He abundance in the presolar material 
(Geiss 1993) and the local ISM (Gloeckler & Geiss 1996). 
However, the abundance in PNe agrees with the predic- 
tions of standard stellar evolution models for stars of mass 
1-1.5 M Q (see the review by Rood et al. 1998). The main 
question is then: if low-mass stars are net producers of 
3 He and return it to the ISM during the PN phase, why 
don't we observe a much higher abundance in Hn regions 
and in the solar system material, as all standard Galactic 
evolutionary models predict (see e.g. Tosi 1996)? 

An interesting solution to this problem involves the 
existence of a nonstandard mixing mechanism (or Cool 
Bottom Processing, hereafter CBP) which operates dur- 
ing the red giant phase of stars with M < 2 Mq. In 
addition to decreasing the amount of 3 He in the stel- 
lar envelope, this process affects the abundances of other 
important elements, including carbon, as first suggested 
by Hogan (1995), Charbonnel (1995) and Wasserburg et 
al. (1995). In particular, the ratio of 12 C/ 13 C in the en- 
velope is predicted to be much lower than in the stan- 
dard case. For a 1 Mq star, the predicted ratio is about 
5 against the standard value of 25-30 in the red gi- 
ant branch (RGB) and of 20-40 in the asymptotic gi- 
ant branch (AGB). However, the discrepancy becomes 
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larger for more massive stars where the 12 C/ 13 C ratio can 
reach ~100 in the AGB phase (Charbonnel 1995; Weiss et 
al. 1996; Forestini & Charbonnel 1997, hereafter FC; van 
den Hoek & Grocnewegen 1997, herafter HG; Marigo 1998; 
Boothroyd & Sackmann 1999, hereafter BS). 

From an observational viewpoint, it is important to 
obtain accurate measurements of the isotopic ratio in 
those PNe where the 3 He abundance has been determined. 
Should these objects show a high value of 12 C/ 13 C, then 
no modifications to the standard stellar models would 
be required. Otherwise, one has to invoke another selec- 
tive process (mixing, diffusion etc.) that operates on some 
isotopes but not on 3 Hc. However, the number of PNe 
with 3 He measurements is small (e.g. Balser et al. 1999), 
whereas the suggested physical processes should be quite 
general and should affect the nucleosynthetic yields of all 
stars of mass less than ~ 2 M Q . Hence, it is critical to 
measure the carbon isotopic ratio in a sample of PNe as 
large as possible. 

The molecular envelopes of PNe have been studied ex- 
tensively at near infrared and millimeter wavelengths (sec 
e.g. Kastncr et al. 1996, Huggins et al. 1996, Bachillcr 
et al. 1997). These observations have shown that massive 
envelopes (>1CP 2 M Q ) containing a rich variety of molec- 
ular species are commonly found around PNe. CO is the 
most widely observed species, and the 12 C/ 13 C isotopic 
ratio has been measured toward several PNe (Bachiller et 
al. 1989, 1997; Cox et al. 1992). These initial studies have 
shown that the 12 C/ 13 C rat io is in the range 10-20. 

Our project consists of two parts. In the first one, we 
have carried out high quality observations of 12 CO and 
13 CO in six PNe that have been searched for 3 He emission. 
In the second run, a larger sample of nebulae with strong 
12 CO line emission has been observed in the 13 CO lines in 
order to determine the isotopic ratio in PNe without 3 He 
measurements. Galli et al. (1997, hereafter GSTP) have 
argued that extra-mixing processes must be at work in 
more than 70% of low- mass stars (M < 2 Mq) in order 
to reconcile the predictions of the Galactic evolution of 
3 He with the observational constraints. We set out this 
experiment to determine the isotopic ratio in a relatively 
large sample of PNe. 

2. Observations 

The observations were carried out with the IRAM 30-m 
telescope at Pico Veleta (near Granada, Spain) during two 
observing runs in November 1996 and May 1997. In the 
first one we observed the six PNe studied by Balser et 
al. (1997). The observations were made simultaneously in 
the J = 2-1 and 1-0 lines of 12 CO. The strongest emitters 
were then observed in the 13 CO lines. In the second run 
we searched for 13 CO J = 2-1 and 1-0 emission in 22 ob- 
jects previously detected in 12 CO by Huggins et al. (1996). 
At the rest frequency of the J = 2-1 and 1-0 12 CO lines 
(near 115 and 230 GHz, respectively) the telescope beam- 



size and efficiency are 24" and 0.6 (115 GHz) and 12" 
and 0.45 (230 GHz). The smaller beam size at 230 GHz, 
together with the small scale structure usually found in 
PNe and with the high intrinsic CO (2-l)/(l-0) line ra- 
tio, lead to observed 2-1/1-0 line ratios in the range 2-5 
(Bachiller et al. 1993), meaning that the J = 2-1 line is 
much more effective for CO searches. The spectrometers 
used were two filter banks of 512 x 1 MHz channels, pro- 
viding a spectral resolution of 1.3 km s _1 in the 12 CO 2-1 
line, and 2.6 km s _1 in the 12 CO 1-0 line. Typical system 
temperatures were around 600 K at 115 GHz and 800- 
1000 K at 230 GHz. The spectra were calibrated with the 
standard chopper wheel method and are reported here in 
units of main-beam brightness temperature (Tmb)- 

3. Measuring the carbon isotopic ratio from 
millimeter wave observations 

Molecular line observations at mm-wavelengths provide a 
powerful method to estimate the 12 C/ 13 C ratio in PNe. In 
particular, the 12 CO/ 13 CO ratio should faithfully reflect 
the atomic 12 C/ 13 C ratio, since the mechanisms which 
could alter the 12 CO/ 13 CO ratio are not expected to be 
at work in PNe. In fact, the kinetic temperature in PN 
envelopes (25-50 K) is high enough that isotopic fraction- 
ation should not operate. Also, selective photodissociation 
is expected to be compensated by the isotope exchange 
reaction 12 CO+ 13 C+ -^ 13 CO+C+ which is faster than 
the 13 CO photodestruction in PN envelopes (e.g. Likkel 
et al. 1988). However, although the J = 1-0 and J = 2-1 
lines of 12 CO have been extensively observed in PNe (e.g. 
Huggins et al. 1996), very few observations of the 13 CO 
lines are available and the value of the isotopic ratio is 
presently poorly known. 

In order to estimate the 12 CO/ 13 CO isotopic ratio, one 
needs to make a number of approximations. First, we as- 
sume that the emitting regions fill the antenna beams in 
the lines of both molecules, or that the filling factor is 
the same (in the case of an extended clumpy medium). 
Second, we assume that the rotational levels are thermal- 
ized at a representative uniform temperature of 25 K (see 
e.g. Bachiller et al. 1997). Thermalization is indeed a rea- 
sonable assumption for 12 CO and 13 CO, since the dipole 
moment is quite small (about 0.1 Debye). Third, if we 
assume that the emission is optically thin for both the 
12 CO and 13 CO lines, then the 12 CO/ 13 CO column den- 
sity ratio is given by the ratio of the integrated intensities. 
We will discuss below the uncertainties introduced by this 
approach in the derived isotopic ratios. 

4. Results 

4-1. PNe with 3 He measurements 

The results of the observations are given in Table 1 where 
we list the PN name, the intensities I21 of the 12 CO and 
13 CO J = 2-1 lines and the derived isotopic ratio. The 
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Fig. 1. CO emission in PNe. Two examples are shown from the sample of PNe with (left) and without (right) 3 He 
measurements. The offset in arcsec from the central PN is also indicated. The stronger J = 2-1 transitions are shown 
in the upper panels, while the J = 1-0 transitions in the bottom panels. The isotopic lines are detected in both cases. 



sample of the six PNe studied in 3 He shows little emission 
in CO: with the exception of NGC 6720, the values of I21 
represent upper limits to the intensity and have been es- 
timated from the line widths deduced from the expansion 
velocities listed in the catalog of Acker et al. (1992). 

Bachiller et al. (1989) detected an extended molecular 
envelope in NGC 6720, the Ring nebula. The CO emission 
reveals a clumpy ring, resembling to that of the ionized 
nebula. We observed the J — 2-1 and 1-0 lines of 13 CO in 
three positions and detected emission in all cases. Fig. 1 
shows the spectra of the four CO and 13 CO lines observed 
toward the strongest peak in the molecular envelope (off- 
set [-40", -20"] from the central position). The measured 
isotopic ratio toward this position, where the S/N is the 
highest, is 12 C/ 13 C = 22. This value is in agreement with 
a previous estimate of Bachiller et al. (1997) and with the 
values obtained toward the two positions with lower S/N 
spectra. Thus, our data provide no evidence for variations 
of the isotopic ratio across the nebula. 

We have also detected a line around the 12 CO J = 1 
frequency in the central position of NGC 6543, but we 



Table 1. CO results for PNe observed in 3 He 



PN Name 


(K km s" 1 ) 


/ 2 i( ld CO) "C/"C 
(K km s" 1 ) 


IC 289 


<0.90 




NGC 3242 


<0.16 




NGC 6543 


<0.19 




NGC 6720 


20.0 


0.9 22 


NGC 7009 


<0.12 




NGC 7662 


<0.12 





failed to detect the J = 2-1 line at relatively low levels. 
This indicates that the line near the J = 1-0 frequency 
is probably not due to CO. It is interesting to recall that 
the H38a recombination line is only separated by 3 MHz 
(7.8 km s" 1 ) from the 12 CO J = 1-0 line. As discussed 
in Bachiller et al. (1992), the H38a line can dominate the 
emission around the 12 CO J = 1-0 frequency in some 
nebulae with little or no molecular gas. We believe that 
this is the case in the central position of NGC 6543. 
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Table 2. CO results for PNe not observed in 3 He 



T mb (K) 




-150 -100 
LSR Velocity (km/s) 



Fig. 2. 13 CO J = 2-1 and 1-0 spectra observed toward 
the central position of the PN identified in the upper left 
corner of each panel. 



PN Name 


/ 2 i( lz CO) 


/ 2 i( ld CO) 




(Kkms- 1 ) (Kkms- 1 ) 




NGC 7027 


448.7 


17.3 


25 T 


NGC 2346 


14.2 


0.58 


23 


NGC 7293 


17.5 


1.8 


9 


NGC 6781 


28.4 


1.72 


20 


M 1-16 


109.0 


32.3 


3 f 


M 1-17 


66.2 


3.2 


22 


M 2-9 


3.3 


2.0 


2 T 


M 2-51 


33.1 


2.3 


15 


M 4-9 


32.2 


1.8 


18 


CRL 2688 


312.5 


112 


3+ 


CRL 618 


182.2 


41.2 


4t 


OH 09+1.3 


5.0 


2.5 


2t 


IC 5117 


19.6 


1.4 


14 


NGC 7008 


1.8 


<0.15 


>12 


NGC 6853 


6.9 


<0.15 


>46 


M 1-13 


22.2 


<1.5 


>15 


BD+30°3639 


4.7 


<1.1 


>4 



t lower limit 



4-2. PNe without 3 He measurements 

The parameters of the detected sources and the derived 
isotopic ratios are listed in Table 2. We have also in- 
cluded in this table data for NGC 2346, NGC 7293, 
NGC 6781, M 4-9, and CRL 618 from previous obser- 
vations by Bachiller ct al. (1989, 1997). The spectra of 
the remaining PNe are shown in Fig. 1 (for M 1-16) and 
in Fig. 2. In four cases (NGC 7008, NGC 6853, M 1- 
13, BD+30°3639) we obtained tentative detections in the 
13 CO lines, or only crude upper limits could be derived. 

The assumption of optically thin emission used to 
determine the isotopic ratio is likely to be accurate in 
the case of evolved nebulae like the Ring, the Helix, 
NGC 2346, etc. In fact, Large- Velocity-Gradient (LVG) 
models confirm that the 12 CO and 13 CO line emission 
is optically thin in these cases. On the other hand, the 
same assumption is not appropriate for the CO lines in 
young objects such as CRL 2688, CRL 618, NGC 7027, 
and M 1-16. In such cases, the derived CO column densi- 
ties and 12 CO/ 13 CO column density ratios represent only 
approximate lower limits. 

One could also have weak emission arising from small 
optically thick clumps very diluted within the 12 CO and 
13 CO beams. This could be the case in some compact CO 
envelopes such as the Butterfly nebula, M 2-9. The CO in 
M 2-9 is concentrated in an expanding clumpy ring which 
has a mean diameter of 6 arcsec (Zweigle et al. 1997). 
Individual clumps in this ring have sizes < 4 arcsec. The 
weakness of the 12 CO and 13 CO lines we observe could be 
due to the important dilution of such small clumps within 
the 30-m beam. The clumps could be optically thick in 
CO, and the reported 12 CO/ 13 CO intensity ratio would 
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just represent a lower limit to the abundance ratio in these 
compact PNe. 

In the case of extended PNe (e.g. NGC 6720, 
NGC 2346, NGC 7293), the assumption of optically thin 
CO 2-1 emission is based on detailed studies of individ- 
ual nebulae which show that the intrinsic 2-1/1-0 line 
intensity ratio is typically >2 (Bachillcr ct al. 1989, 1993; 
Huggins et al. 1996). Such high ratios imply that the emis- 
sion is optically thin, and that the excitation temperature 
is Tex^lO K. In this approximation, and assuming homo- 
geneous excitation conditions along the line of sight, the 
column density in the J=2 level is proportional to the 
observed 2-1 line intensity. Moreover, the relative popu- 
lation of the J=2 level is quite insensitive to the value of 
the excitation temperature for the typical conditions of 
PNe. The total CO column densities determined in this 
way are correct for T ex in the range from 7 to 77 K, and 
are within a factor of 2 over the range from 5 to 150 K, 
which should cover most PNe. In any event, because of 
the similar dipole moments, the excitation temperature 
should be the same for both CO and 13 CO. Then, if the 
excitation conditions do not vary along the line of sight, 
the isotopic ratios are independent of T ex . 

The major source of uncertainty for these extended 
nebulae is thus related to the observational procedures. 
The calibration of our observations is accurate within 20%, 
but small differences in the filling factor for the differ- 
ent spectral lines could increase the uncertainty of the 
measured isotopic ratios up to a factor ~1.5. As Table 2 
shows, in the extended PNe where 13 CO is well detected 
(NGC 6720, NGC 2346, NGC 7293, NGC 6781, M 4-9, 
and M 2-51) the values of the isotopic ratios are in the 
range 9 to 23, to within a factor of ~ 1.5. Such values are 
thus significantly lower than the Solar System value of 89, 
and appear to be consistent with those derived in AGB 
stars, as we shall discuss in Sect. 6.1. 

5. Mass estimate of the progenitors of the PNe 

In order to compare the observed isotopic ratios with the 
predictions of stellar evolutionary models, we must esti- 
mate the mass of the progenitors stars of the PNe. For the 
present samples, we have followed the procedure adopted 
in GSTP, consisting of the following steps. 

(1) We assign to each nebula the best measured dis- 
tance (or an average of available distances), the H/3 flux, 
the Hen flux at 4686 A, the angular size, and the B and V 
stellar magnitudes. Using these parameters, we calculate 
the effective temperature and luminosity via the Zanstra 
method (Kaler 1983). 

(2) By locating the central star in the log T e ff-log L 
plane, we derive its mass (Mqs) from comparison with a 
set of evolutionary tracks (Stanghcllini & Renzini 1993). 

(3) Using the initial mass-final mass relation, we com- 
pute the progenitor mass, that is the stellar mass on the 
main sequence (Mms)- 
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The stellar properties adopted for the PNe detected 
in 13 CO and the derived values of the progenitor mass 
are given in Table 3. Details on the individual objects are 
given in the Appendix. 

Let us examine the uncertainty involved in the final 
mass calculations. Estimates of the stellar temperature 
and luminosity given in Table 3 are affected by errors 
in magnitudes, fluxes, diameters and extinctions. How- 
ever, these quantities are usually determined with good 
accuracy (~ 5 — 20%), so that the uncertainty in the de- 
rived mass of the central stars does not exceed ~ 15%, or 
~ 0.02 Mq. The values given in the table do not include 
errors on the distances to the PN, which can be intrinsi- 
cally high (up to 50%) but are difficult to estimate on an 
individual basis. 

To infer the main sequence masses, we have used the 
initial mass-final mass relation given by Hervig (1996). 
This relation differs from that of Weidemann (1987) 
adopted in GSTP. We preferred Hervig's prescription since 
it is based on reliable observations of cluster white dwarfs, 
although the formal errors on the final mass arc still sub- 
stantial, and can amount to about 0.1 Mq. Since we de- 
rive initial masses from final masses, the errors on the 
former quantity can be even larger. Quantitatively, we as- 
sign a formal error to the main sequence mass of AMms — 
1.5 M Q for low values of the initial mass (Mms < 2 Mq), 
and a smaller error (AMms — 0.75 Mq) for higher masses. 
This difference is due to the change of the slope of the ini- 
tial mass - final mass relation at about 2 M©: smaller 
masses are more sensitive to the adopted relation, and the 
uncertainty is correspondingly larger. 

6. Implications 

6.1. Implications on stellar nucleosynthesis 

How do we interpret our results on the 12 C/ 13 C ratios 
in the framework of stellar nucleosynthesis? To help an- 
swering this question, we combine the information pro- 
vided by the observed 12 C/ 13 C isotopic ratios with the 
mass estimates of the progenitors of the PNe, and with the 
predictions of some representative stellar nucleosynthesis 
models. 

Since the formation of a PN takes place at the end of 
the AGB phase, the significant comparison is between the 
observed abundances and those predicted for the stellar 
ejecta at the AGB tip. Unfortunately, no stellar nucleosyn- 
thesis models up to these late phases are available in the 
literature for stars experiencing deep-mixing. For example, 
the calculations of Boothroyd & Sackmann (1999) allow 
to derive the mass dependence of the 12 C/ 13 C ratio on 
the stellar surface at the tip of the red giant branch both 
for the standard case and in the presence of extra-mixing. 
The standard models indicate that the isotopic ratio has 
approximately a constant value of 20-23 between M = 2 
and 4 Mq, and then increases steadily at lower masses 
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Table 3. Properties of individual PNe detected in 13 C0 



# 


PN Name 


Log T efI 


Log L 


D 


M C s 


M MS 






(K) 


(Le) 


(kpc) 


(M & ) 


(M e ) 


1 


NGC 7027 


5.28 ±0.02 


3.66 ±0.03 


0.88 


0.67 ±0.03 


2.6 


2 


NGC 2346 


4.76 ±0.01 


3.80 ± 0.04 


0.90 


0.63 ±0.02 


2.3 


3 


NGC 7293 


5.03 ±0.01 


2.03 ±0.03 


0.28 


0.60 ±0.03 


2.6 


4 


NGC 6781 


4.99 ± 0.03 


2.94 ±0.09 


1.57 


0.55 ±0.02 


1.0 


5 


M 1-16 


4.91 ±0.03 


3.36 ±0.14 


5.43 


0.56 ±0.02 


1.0 


6 


M 1-17 


4.97 ±0.03 


2.89 ±0.11 


7.36 


0.55 ±0.05 


1.0 


7 


M 2-9 


> 4.65 


> 2.54 


1.32 


C) 


> 2 


8 


M 2-51 


5.06 ±0.06 


2.20 ±0.15 


1.92 


0.63 ± 0.09 


2.3 


9 


M 4-9 


> 4.85 


> 0.18 


1.71 


(*) 




10 


CRL 2688 










> 1.3 


11 


CRL 618 


> 4.48 


> 1.11 


2.81 


(*) 




12 


OH 09+1.3 






8.0 






13 


IC 5117 


4.95 ±0.03 


3.31 ±0.12 


2.10 


0.56 ±0.02 


1.0 



(*) no Hen, no Zanstra analysis possible 




1 2 3 4 5 

Progenitor Mass (M ) 



Fig. 3. 12 C/ 13 C versus progenitor mass of the PNe de- 
tected in CO. The filled triangles are the measured PNe, 
labelled as in Table 3. The uncertainty on the mass 
is described in the text. The error bar on the derived 
12 C/ 13 C ratio is shown in the lower right corner. The three 
curves are the predicted 12 C/ 13 C ratios in the ejecta of 
stars at the tip of the asymptotic giant branch (AGB). 
The dotted curve is from van den Hoek & Grocnewegen 
(1997), the dashed curve is from Forestini & Charbonnel 
(1997), and the long-dashed curve is from Marigo (1998). 
All curves are for solar metallicity. 

up to about 28-30 with a small dependence on the stellar 
metallicity. In the case of extra mixing, the 12 C/ 13 C ratio 
displays a sharp drop below ~ 2 M©, reaching values of 



5-10 at ~1 M©. Similar results have also been obtained 
by Charbonnel (1994) and Denissenkov & Weiss (1996). 

In Fig. 3 we show the distribution of the measured 
12 C/ 13 C ratios in 11 PNe, together with the theoretical 
values in the ejecta of stars at the tip of the AGB phase. 
The different curves refer to the models computed in the 
standard case with no mixing and solar metallicity by FC, 
HG, and Marigo (1998). Fig. 3 indicates that over the mass 
interval between 1.5 and 4 M Q the predictions of the HG 
and Marigo models are in good agreement with each other 
and indicate a roughly constant value of 12 C/ 13 C ^100. 
The models of FC show significantly lower ratios at 3 and 
4 Mq, but have the same qualitative behavior. These re- 
sults are well understood in terms of the nucleosynthesis 
occurring during the thermally pulsating AGB phase. In 
fact, a major phase of 12 C enrichment of the convective 
envelope results from the penetration of the convective 
tongue during thermal pulses (see the discussion in, e.g., 
FC and HG). At the same time, 13 C is partially burnt 
through the 13 C(a,n) reaction at the bottom of the inter- 
shell region during the inter-pulse phase, as first suggested 
by Straniero et al. (1995). The combination of these two 
effects accounts for the high 12 C/ 13 C ratio in this mass 
range. On the other hand, more massive AGB stars ex- 
perience hot bottom burning that progressively leads the 
12 C/ 13 C ratio close to its equilibrium value (~4-5). 

Despite the uncertainties in our measurements, it ap- 
pears that the standard models produce 12 C/ 13 C in excess 
of the values observed in our objects. The only exception 
is for the two lowest mass stars (0.9 and 1.0 M Q ) com- 
puted by HG which agree with the ratios measured in 
Ml-17 and NGC 6781. Although Fig. 3 indicates a marked 
discrepancy between observed and theoretical values, one 
should be aware of the sensitivity of the predicted yields 
on the model assumptions. In particular, the isotopic ra- 
tios (not only 12 C/ 13 C ) depend rather strongly on the 
adopted mass loss rate during the AGB phase, on the third 
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dredgc-up efficiency, and on stellar metallicity. Their com- 
bined effects have been thoroughly discussed by FC and 
HG and need not to be repeated here. In general, a shorter 
phase of mass loss and/or lower mass loss rates tend to 
decrease the surface 12 C/ 13 C ratio. However, a more effi- 
cient convective penetration results in a higher pollution of 
the envelope from pulse to pulse and a higher 12 C/ 13 C ra- 
tio. Finally, stars with initial metallicities lower than solar 
yield higher 12 C/ 13 C ratios. 

Returning to Fig. 3, we see that several PNe lie in the 
region of the diagram with low 12 C/ 13 C and low mass. 
For these objects we should expect the standard models 
to become inaccurate insofar as they neglect the possible 
occurrence of mixing mechanisms of non-convective origin 
which can alter the composition of the stellar ejecta. Low 
values of the carbon isotopic ratio have also been mea- 
sured in field population II stars and in globular cluster 
giant (e.g. Sncdcn et al. 1986; Gilroy & Brown 1991; Pi- 
lachowski et al. 1997) and have provided the motivation 
to introduce extra-mixing processes in the standard evo- 
lution (e.g. Charbonnel 1995). Recently, Charbonnel & do 
Nascimento (1998) find that more than 90% of a sample 
of 191 field and cluster red giants presents carbon isotopic 
ratios inconsistent with those predicted by standard nu- 
cleosynthesis. 

The observational situation in AGB stars is less clear, 
because of the extreme sensitivity of the determination of 
isotopic ratios on excitation temperatures and model at- 
mospheres. For example, Greaves & Holland (1997) have 
measured a 12 C/ 13 C ratio varying between 12 and 36 in a 
sample of 9 carbon stars with high mass-loss rates (~ 10~ 5 
Mq yr _1 ), in accordance with similar results previously 
obtained by Wannier & Sahai (1987) on seven other car- 
bon stars. On the other hand, Lambert et al. (1986) mea- 
sured the 12 C/ 13 C ratio in 30 cool carbon stars and found 
values between 30 and 70. More recently, Ohnaka & Tsuji 
(1996) report ratios between 20 and 30 for 24 N-type car- 
bon stars in common with the sample of Lambert et al. 
(1986). These values are about a factor of 2 smaller than 
those derived by Lambert et al. (see however the discus- 
sion in de Laverny & Gustafsson 1998 and Ohnaka & Tsuji 
1998), but are consistent with those found for our PN sam- 
ple. If these samples are representative of the population 
of both carbon stars and PNe, the low 12 C/ 13 C ratios 
provide strong indication that most stars should experi- 
ence some extra-mixing and deplete 12 C with respect to 
13 C. According to this conjecture, it appears that in order 
to match the observations, more realistic stellar models 
should include additional physical processes to avoid too 
large 12 C/ 13 C ratios throughout the latest stages of stel- 
lar evolution. 

6.2. Implications on Galactic chemical evolution 

Theoretical models predict that stars with low values of 
12 C/ 13 C should also have low 3 Hc abundances. In GSTP, 
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Fig. 4. Evolution of 12 C/ 13 C and 3 He as a function of time 
for the solar neighborhood and as a function of Galac- 
tocentric radius at the present time. Standard case with 
no extra mixing. The three curves show the results ob- 
tained using the yields for intermediate-mass stars by FC 
(solid), by HG (dotted), and by BS (dashed). The obser- 
vational data for 12 C/ 13 C are from: Anders & Grevesse 
(1989); Henkel et al. (1980); Gardner & Whiteoak (1981); 
Henkel et al. (1982); Henkel et al. (1985), Langer & Pen- 
zias (1993); Wouterloot & Brand (1996). The data on 3 He 
are from Geiss (1993), Gloeckler & Geiss (1996) and Rood 
et al. (1995). 

we estimated that to obtain consistency between the ob- 
served 3 He abundances and chemical evolution models, 
the majority of the Galactic PNe (>70%) should have a 
12 C/ 13 C ratio well below the standard value. In this sec- 
tion we test this suggestion using the constraint provided 
by the C/ 13 C observations in PNe. Chemical evolution 
models offer the adequate tool to follow simultaneously 
the evolution of 12 C/ 13 C and 3 He over the Galactic life- 
time. 

In Figs. 4 and 5 we show the evolution of 12 C/ 13 C and 
3 He as a function of time in the solar neighborhood 
(Galactocentric radius R = 8 kpc), and as a function 
of R at the present time (t = 13 Gyr), as predicted by 
models described by Dearborn et al. (1996) and Sandrclli 
et al. (1998). The two figures correspond to different as- 
sumptions on the fraction of low-mass stars experiencing 
deep mixing in the RGB phase. In both cases, the Galactic 
evolution of the various isotopes has been computed as- 
suming metallicity-dependent yields. In the upper panels, 
the solid line corresponds to models adopting the 12 C and 
13 C yields from BS for stars in the mass range 0.8 - 2.5 
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Fig. 5. Evolution of 12 C/ 13 C and 3 He as a function of time 
for the solar neighborhood and as function of Galactocen- 
tric radius at the present time. Case with mixing in 90% 
of low mass stars. Symbols and observational data are as 
in Fig. 4. 

Mq, from FC in the range 3-6 M©, and from Woosley & 
Weaver (1995) for massive stars. The dotted and dashed 
lines display the results obtained using the same yields as 
before, except for stars in the mass range 3-6 M©, where 
the HG and BS yields, respectively, are adopted. The cal- 
culations by BS are the only ones that provide the 12 C 
and 13 C yields for both the standard and the CBP cases. 
For 3 He, the standard yields are taken from Dearborn ct 
al. (1996) and the CBP yields from BS. 

The behaviour of 3 He is typical of the elements pro- 
duced mainly by low-mass, long-lived stars. 12 C is a pri- 
mary element produced by stars of any mass and avail- 
able for Galactic enrichment already after the explosions 
of the first massive objects. 13 C is mostly secondary and its 
bulk abundance is due to intermediate-mass stars. There- 
fore, the enrichment of this element in the ISM occurs 
later than that of 12 C. This delay causes the decrease of 
12 C/ 13 C visible in Figs. 4 and 5. The positive radial gra- 
dient of 12 C/ 13 C is also due to the different mass intervals 
of stars that are producers of the carbon isotopes, and to 
the higher star formation activity in the inner Galactic 
regions, independent of the stellar initial mass function. 
All these effects combined make the relative proportion of 
13 C over 12 C producers increasingly higher for decreasing 
Galactocentric distances. 

Fig. 4 shows the standard case with no deep mixing, 
whereas the results with CBP in 90% of stars randomly 
chosen in the mass range M < 2.5 M© are displayed in 



Fig. 5. The hatched regions represent the range of abun- 
dances resulting from 50 different cases for the random 
sorting. As already discussed by GSTP, the overall en- 
richment of 3 He is very sensitive to the assumptions on 
deep mixing. On the contrary, the variation of 12 C/ 13 C is 
marginal, since most of 12 C and 13 C is produced by stars 
more massive than 2.5 M© which are not affected by deep 
mixing processes. Thus, altering the amount of 12 C and 
13 C released by low- mass stars leads to smaller variations 
in the resulting abundances. For these reasons, the Galac- 
tic evolution of 3 He is an excellent indicator of the fraction 
of low-mass stars which should undergo CBP, but that of 
12 C/ 13 C is not. From Figs. 4 and 5, one can see that the 
fraction required to reproduce the observed pattern of 3 He 
abundances is up to 90%, in agreement with the finding 
of Charbonnel & do Nascimento (1998). 

As for 12 C/ 13 C, the predicted abundance distributions 
with time and with Galactocentric distance are roughly 
consistent with the data regardless of the occurence of 
extra-mixing. The 12 C/ 13 C ratio in the Sun cannot be 
accurately reproduced, and the radial gradient is flatter 
than the observed one (as deduced from measurements in 
molecular clouds), independently of the adopted percent- 
age of stars undergoing CBP. Using the yields of FC and 
BS leads to a better fit of the solar value and to the present 
distribution of 12 C/ 13 C in the inner Galactic regions. On 
the other hand, the models of HG reproduce the local ISM 
ratio. In all cases, the fraction of low-mass stars experi- 
encing deep mixing does not affect significantly the overall 
results of the chemical evolution models. Of course, deep 
mixing alters the values 12 C/ 13 C in low-mass stars, but 
the Galactic evolution of 12 C and 13 C is mainly governed 
by stars in which this process is not expected to take place. 

7. Prospects for future observations 

One of the main goals of our study was the determina- 
tion of the isotopic ratio in NGC 3242, the best object for 
testing the mixing hypothesis since the 3 He measurements 
are very reliable and the progenitor mass is in the range 
where extra mixing should occur (see also the discussion 
in Rood et al. 1998). However, the high sensitivity of our 
observations and the fact that we searched for CO emis- 
sion across the whole circumstellar shell imply that CO 
is absent in the nebula. We therefore conclude that the 
isotopic ratio cannot be determined in this important ob- 
ject by means of millimeter line observations. The lack of 
a molecular envelope around NGC 3242 can be a conse- 
quence of the low-mass of the progenitor. It is well known 
that the transition from the AGB phase to the PN stage 
is slow for low-mass stars and therefore the gas remains 
exposed to ionizing and dissociating radiation for a longer 
time than in the case of more massive stars. 

This behavior is consistent with the fact that the most 
massive object of our sample, NGC 6720, indeed has a rich 
molecular envelope. Unfortunately, NGC 6720 has an esti- 
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mated progenitor mass of <~ 2 M© , quite close to the limit 
where the nonstandard mixing mechanism is expected to 
significantly affect the isotopic ratio. Thus, the derived 
ratio is consistent with both standard and nonstandard 
evolutionary models. The lesson is clear: there is a trade- 
off between the need to select low-mass PNe in order to 
discriminate between theoretical models and the higher 
detection rate of molecular envelopes around more mas- 
sive objects. The determination of isotopic ratios may be 
a tricky business, but future observations should take up 
such a challenge. 

The method of using mm-wave transitions is not the 
only one adequate for measuring the isotopic abundance 
in PNe. Clegg (1985) first pointed out the possibility of 
using the CIII] multiplet near 1908 A to obtain a di- 
rect estimate of the isotopic ratio in the ionized gas of 
PNe. Recently, Clegg et al. (1997) have successfully de- 
tected the extremely weak isotopic line 13 C \/2^o~\/2^ 
in two PNe (plus a tentative detection in a third one), 
using the HST Goddard High Resolution Spectrograph. 
The derived 12 C/ 13 C ratios arc 15±3 in NGC 3918 and 
21±11 in SMC N2. In either case, no measurement at 
mm-wavelengths has been made to independently check 
the derived values. However, as Clegg et al. point out, the 
two types of measurements are complementary since they 
cannot be performed on the same objects: the UV tran- 
sitions require high excitation conditions, while the op- 
posite is true for the mm-wave lines studied by us. Con- 
sidering the negative results of the CO search, it would 
be extremely important to perform the HST observations 
towards NGC 3242. 



8. Conclusions 

We have observed a sample of 28 PNe in 12 CO and 13 CO. 
We have determined the isotopic ratio in 14 objects, 9 
of which are new detections, and obtained robust upper 
limits in 4 other PNe. The observations reported in this 
paper improve the available data on the 12 C/ 13 C ratio in 
PNe. The results confirm and extend earlier measurements 
which indicate values between 5 and 30, below those pre- 
dicted by current AGB models. According to these cal- 
culations, the value of 12 C/ 13 C in the ejecta of stars in 
the mass range between ~ 1.5 and ~ 4 M© is of the or- 
der of 100. In order to agree with the observations, we 
suggest that stellar nucleosynthesis models should include 
the effects of extra-mixing to avoid the overproduction of 
12 C with respect to 13 C not only during the RGB phase, 
but also in the latest phases of evolution. Our estimates 
of the mass of the progenitor stars support this conclu- 
sion, since all PNe considered here originated from stars 
with mass less than ~ 3 M©. The proposed extra-mixing 
mechanism is expected to be at work in this mass range 
(Hogan 1995, Charbonnel 1995) and represents at present 
the most likely solution to the so-called " 3 He problem" . 



Out of the six PNe with measured 3 He abundances, 
only for one (NGC 6720) it was possible to estimate the 
12 C/ 13 C ratio. The remaining five objects show weak or 
no 12 CO emission. In the case of NGC 3242 where the 
3 He measurements are the most reliable, the absence of 
CO emission in the envelope precludes a direct test of the 
predictions of stellar nucleosynthesis using millimeter line 
observations. We have suggested an alternative method to 
determine the isotopic abundance in this critical object. 

We have used chemical evolution models to follow the 
variation with time and Galactocentric distance of the 3 He 
abundance and 12 C/ 13 C ratio. Whereas the model predic- 
tions for 3 He are inconsistent with the observed abundance 
in the Solar System and in the interstellar medium unless 
extra-mixing is operative in at least 90% of low-mass stars, 
the evolution of 12 C/ 13 C is remarkably insensitive to the 
effects of extra- mixing. The Galactic evolution of 12 C and 
13 C is mainly governed by intermediate and massive stars, 
in which extra-mixing is not expected to occur. 
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Appendix A: Notes on individual PNe 

Here, we provide detailed information on the physical 
properties of the PNe of our sample. Unless otherwise 
noted, B and V magnitudes are from the Acker et 
al. (1992) catalogue. Fluxes, extinction constants, and di- 
ameters are from Cahn et al. (1992, hereafter CKS). 
NGC 7027: The Zanstra analysis yields logT cff = 5.28 ± 
0.02 and log L/L Q = 3.66 ± 0.03. We have used the ex- 
pansion distance from the radio image of Masson (1989). 
The derived mass of the central star is Mcs = 0.67 M©. 
NGC 2346: The distance, D = 0.9 kpc, is derived from 
the average extinction distance (Pottasch 1996). We ob- 
tain logT off = 4.76 ± 0.01 and logi/L© = 3.80 ± 0.04. 
The mass of the central star results Mcs = 0.63 M Q . 
NGC 7293: Many distance measurements are available 
for this nebula. Pottasch (1996) quotes a very uncertain 
(error up to 100%) parallax distance of D = 0.19 kpc. 
A more reliable value of D = 0.28 kpc comes from av- 
eraging several individual distances (Pottasch 1996). We 
obtain logT cff = 5.03 ± 0.01 and logi/L© = 1.80 ± 0.03 
or 2.03 ± 0.03, using the parallax and average distances, 
respectively. The corresponding central star masses are 
Mcs — 0.675 and 0.60 M . The progenitor masses are 
2.6 and 2.0 M©. Given the large error on the parallax, we 
prefer the latter value for the initial mass of NGC 7293. 
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NGC 6781: We obtain logT eff = 4.99 ± 0.03. By using 
the statistical (CKS) and the average of individual (Acker 
et al. 1992) distances, we obtain respectively logi/T = 
2.24 ± 0.09 and 2.94 ± 0.09. In both cases, the mass of 
the central star coincides with the evolutionary track for 
M C s = 0.55 M Q . 

M 1-16: By using D = 5.45 kpc (CKS method with a 
newly measured diameter), wc find logL/ 'Lq = 3.36 ±0.14 
and log T c ff = 4.91 ±0.03. The inferred mass of the central 
star is M cs = 0.56 M Q . 

M 1-17: With a distance of D = 7.36 kpc (from CKS 
method), we obtain logT ff = 4.97 ±0.03 and log L/ Lq = 
2.89 ±0.11. The location on the HR diagram is just below 
the 0.55 M Q . 

M 2-9: We calculate the distance with the method of 
CKS, but we evaluate the effective diameter from the Ha 
image (6 = 34.75 arcsec). The new distance to this nebula 
is D = 1.32 kpc. No Hen flux has been detected. Thus, the 
Zanstra analysis gives lower limits to the temperature and 
luminosity. We obtain log Toff > 4.65 and log L/Lq > 2.54 
from the hydrogen recombination lines. To date, nothing 
has been published on this nebula to set better constraints 
on the mass of the central star. 

M 2-51: The Zanstra analysis gives logT cff = 5.06 ±0.06 
and logT/L Q = 2.20 ± 0.15. The luminosity is estimated 
using a distance D — 1.92 kpc, derived with the CKS 
method and a new estimate of the diameter. The resulting 
final mass is Mcs = 0.63 M . 

M 4-9: The distance evaluated with the newly mea- 
sured angular size based on the Ha image by Schwarzet 
al. (1992) is D = 1.71 kpc. For the central star, only a pho- 
tographic magnitude exists. Thus, the resulting Zanstra 
analysis is quite uncertain. The Hp flux is from Acker et 
al. . The Hen flux has never been measured. We obtain 
log T c ff > 4.85 and logT/L Q > 0.18. A value of the mass 
of the central star cannot be derived. 
CRL 2688: Catalogued as a "possible planetary neb- 
ula" , the famous Egg nebula is actually a proto-PN. The 
standard analysis via Zanstra method is not feasible, as 
the nebula is still very thick to optical radiation. There 
are two mass estimates in the literature, Mms = 1-3 or 
2.7 Mq respectively, which depend on the assumed mass 
loss rate (Sahai et al. 1998). 

CRL 618: The Hell flux is not available. Therefore, we 
cannot derive a Zanstra temperature. Hydrogen recom- 
bination lines give logT c ff > 4.48 and a luminosity of 
\ogL/LQ > 1.11, if D = 2.81 kpc is used. This distance 
has been derived with the CKS method and diameter mea- 
sured by Manchado et al. (1996). No mass determination 
is possible. 

OH 09+1.3: Another proto-PN, with no published mass 
to date. 

IC 5117: The distance is the average of the extinction 
distances quoted in Acker et al. (1992). We calculate 
logTeff = 4.95 ± 0.03 and log L/Lq = 3.31 ± 0.12. The 
derived mass of the central star is Mcs = 0.56 M Q . 



M 1-59: Since no stellar magnitudes are available for this 
object, we calculate the effective temperature with the 
"crossover" method (Kaler 1983). We use the statistical 
distance D — 2.50 kpc from CKS, based on a new diameter 
measured by Manchado et al. (1996) and obtain logT ff = 
5.07 ± 0.01 and logi/T Q = 1.37 ± 0.02. The mass of the 
central star mass is Mcs = 0.85 M Q . 
NGC 6853: The parallax distance quoted in Pot- 
tasch (1996) is of the best quality, so we can confidently 
use this value to evaluate the luminosity. From the Zanstra 
analysis we derive logT e ff = 5.14 ± 0.03 and logT/L© = 
2.52 ± 0.07, and a central mass M C s = 0.67 M 
BD+30°3639: The Hen flux is not detected. Thus, we 
obtain lower limits to the luminosity and temperature: 
logT c ff > 4.68 and logL/ Lq > 2.70, using the hydro- 
gen recombination lines and the average of extinction dis- 
tances from Acker et al. (1992). No mass determination is 
possible. 

NGC 7008: We calculate logT cff = 4.98 ± 0.01 and 
logT/L Q = 3.76 ± 0.13. The distance used is the average 
of the extinction distances quoted in Acker et al. (1992). 
The derived central mass is Mcs = 0.63 M Q . 
M 1-7: The Zanstra analysis gives logT cff = 5.05 ± 0.04 
and logL/ Lq = 2.43 ± 0.15 when using D = 3.47 kpc, 
as derived from the diameter of Manchado et al. (1996) 
and the method of CKS. The mass of the central star is 
M cs = 0.56 M Q . 

M 1-13: Since no stellar magnitudes are available for this 
object, we calculate the effective temperature with the 
"crossover" method (Kaler 1983). We use the statistical 
distance D — 5.32 kpc from CKS to obtain logT e ff = 
5.35 ± 0.03 and logL/ Lq = 3.70 ± 0.26. The central star 
mass is M C s = 0.71 M Q . 
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